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Abstract 

We determine characteristic timescales for the viscous damping of r-mode oscihations in rapidly 
rotating compact stars that contain quark matter. We present results for the color-flavor-locked 
(CFL) phase of dense quark matter, in which the up, down and strange quarks are gapped, as well 
as the normal (ungapped) quark phase. While the ungapped quark phase supports a temperature 
window 10® K < T < 5 x 10^ K where the r-mode is damped even for rapid rotation, the r-mode 
in a rapidly rotating pure CFL star is not damped in the temperature range lO^'^ K< T < 10^^ 
K. Rotating hybrid stars with quark matter cores display an instability window whose width is 
determined by the amount of quark matter present, and they can have large spin frequencies outside 
this window. Except at high temperatures T > 10^^ K, the presence of a quark phase allows for 
larger critical frequencies and smaller spin-periods compared to rotating neutron stars. If low-mass 
X-ray binaries contain a large amount of ungapped or CFL quark matter, then our estimates of 
the r-mode instability suggest that there should be a population of rapidly rotating binaries at 
1/ > 1000 Hz which have not yet been observed. 

PACS numbers: 26.60.-hc, 24.85.-hp, 97.60.Jd 
Keywords: r-mode, viscosity, quark matter, CFL 



* Electronic address: jaikumarOimsc . res . in 

t Electronic address: grupakOu. washington.edu. Permanent address: Department of Physics and Astron- 
omy, Mississippi State University, Mississippi State, MS 39762 
■'■ Electronic address: steineraOpa.msu.edu 

1 



I. INTRODUCTION 



Neutron stars are astrophysical test-beds for the physics of strongly-interacting matter 
in regimes of density that are not presently accessible either in terrestrial experiments or by 
numerical studies of lattice-regularized QCD. Therefore, connecting theoretical advances in 
the knowledge of the equation of state (EoS) of dense matter to astrophysical observations 
of neutron stars is essential. The EoS, which plays an important role in determining a 
compact star's properties and in its evolution with time, is constrained by laboratory data 
only up to densities of about 2pAr, where pn ~ 2.6 x 10^^ g cm~^ is the saturation density 
of nuclear matter. At asymptotically high densities where perturbative calculations in QCD 
are reliable, quark matter is conjectured to be in a color superconducting state commonly 
called the color-flavor-locked (CFL) phase [l|. Meanwhile, models of low-energy QCD at 
moderate baryon chemical potentials /xb ^ 1 GeV in a restricted temperature domain T < 
Tc ~ 50 MeV predict a rich patchwork of color superconducting phases in the QCD phase 
diagram 0, i, B, S, 9, with important consequences for the properties of self-bound quark 
stars or neutron stars whose core density is large enough for deconfined quark matter to 
exist there. We refer to the former as strange stars and the latter as hybrid stars throughout 
this work. 

The effects of color superconductivity on the EoS are small, though potentially important 
for the surface structure of strange stars a i a, or the nuclear-qua;k interface in hybrid 
stars More interesting are its effects on transport phenomena, as pairing among the up. 



down and strange quarks breaks the SU{3) global chiral symmetry spontaneously [l|, leading 
to (pseudo-) Nambu-Goldstone modes that determine the material response to external 



perturbations. See the recent review Ref. for the effective theory of these Nambu- 



Goldstone bosons. The pairing pattern that leaves the maximal residual symmetry in the 
three-flavor case is the CFL phase. In this phase, the up, down and strange quarks of 
all three colors participate in Cooper pairing and the order parameter reduces the QCD 
symmetry SU{3)l x SU{3)ji x SU{3)c to the diagonal subgroup SU{3)l+r+c locking the 
color and flavor symmetries. Apart from the 8 pseudo-Nambu-Goldstone bosons associated 
with the spontaneous breaking of chiral symmetry, other light modes in this phase are 
the exactly massless Nambu-Goldstone boson (the superfiuid phonon), which is a long- 
wavelength fluctuation of baryon number, and the massless photon, which is effectively a 
combination of the vacuum photon and a gluon. 

Studies of transport coefficients reveal that the thermal conductivity and specific heat of 
the CFL phase at temperatures much smaller than Tc, the critical temperature for pairing, is 
controlled by the phonon and the photon, while the electrical conductivity is determined by 
thermally excited electron-positron pairs [12]. In the CFL phase the long-term cooling of the 
star is controlled by neutrino emission from the scattering and weak decay of phonons [l3| . 
These novel transport phenomena imply that neutrino transport and neutron star cooling 
curves have the potential to distinguish between the numerous possible patterns of Cooper 



pairing that might occur at stellar densities [IJ, ll5|, ll6|, Il7 



Another set of transport coefficients that are modified by Cooper pairing are the bulk 
and shear viscosities. Bulk viscosity can lead to energy dissipation when the rate of chemical 
equilibration in matter is comparable to the driving frequency of the volume perturbation. 
Shear viscosity determines the relaxation of momentum components perpendicular to the 
direction of fluid flow. Together, they characterize the material response to compressional 
and shearing forces that are typically present in rotating and pulsating compact stars, and 
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serve to damp out large-amplitude oscillatory modes of the fluid. Moreover, since the fluid 
inside a star is self-gravitating, these oscillations can couple to metric perturbations, causing 
the star to emit gravitational waves. Indeed, the primary physical interest in exploring 
viscosities of color superconducting quark matter is to assess their role in damping out 
the so-called r-mode, a particular oscillatory mode of the fluid that dissipates the star's 
rotational energy by efficiently coupling the angular momentum of the star to gravitational 
waves [n, [m, . The damping timescale, which depends sensitively on the viscosity, and 
hence on the low-energy modes of the color superconductor determines if the r-mode will 
be driven unstable. The unsuppressed growth of the r-mode eventually saturates due to 
non-linear effects 2]J since energy can be transferred from r-modes to "inertial" oscillation 
modes which do not couple to gravitational waves. However, the saturation amplitude can 
still be large enough to emit gravitational wave signals that can be detected directly by 
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The observational context of r-modes is that neutron stars do not typically spin at rates 
near the maximum allowed frequency, the Kepler frequency Qk- r-modes offer a possi- 
ble explanation of this fact: in rotating neutron stars, these modes lose energy through 
gravitational waves, which carry away angular momentum from the star and act as braking 
radiation. We mention two distinct settings in which r-modes can be instrumental in slowing 
the neutron star's rotation speed. 

(a) In the newly-formed isolated neutron star in the aftermath of a core-collapse su- 
pernova: Whether r-modes are relevant for neutron stars at birth is difficult to determine 
because the number of rapidly spinning pulsars that can be detected close to their birth is 
a negligible fraction of the total detectable sample. Limited by this selection effect, most 
inferences of rotation rates of neutron stars rely on a model of time evolution of radio pul- 
sars from birth (neglecting any contribution from r-modes), with conflicting results on the 



birth spin-periods, ranging from milliseconds [23| to hundreds of milliseconds [2J]. An in- 
dependent method of determining the initial spin period based on correlating a supernova's 
observed X-ray luminosity (or upper limits thereof) with the rotational energy loss of the 
neutron star suggests that a population of millisecond pulsars at birth is ruled out (25| . 
On the other hand, simulations of Type II supernovae including rotational effects predict 
relatively fast initial spin periods [26|, which is difficult to reconcile with observed young 
and slowly rotating pulsars (like the Vela and Crab pulsars). Thus, r-modes may play an 
important role in neutron stars at birth, and the impact of r-modes depends critically on the 
composition of the neutron star core. The critical temperature for color superconductivity 
being of the order of a proto-neutron star's temperature, it is important to understand the 
r-mode in the quark superfluid as it can determine the initial spin-down evolution of the 
strange star or hybrid star. 

(b) In old neutron stars spun-up by accretion in binary systems: Millisecond pulsars 
(typically, pulsars with spin-periods ~ 10 ms) are old and most of them are observed to be 
in low-mass X-ray binary systems (LMXBs) |27|]. The fast rotation is thought to result from 
spin-up due to steady accretion from a low-mass companion star. This picture has recently 
received conflrmation through the observation of pulsations in LMXBs. Spin frequencies 
in LMXBs, as determined either from X-ray burst oscillations or from kHz quasi-periodic 
oscillations, are typically between 300-640 Hz, well below the maximum rotational rate (see 
Ref. [Isl for a recent review). The upper limit on the rotation rate of LMXBs has been 



recently revised by the detections of objects spinning at 716 Hz [29| and 1122 Hz |30|. How- 



ever, the statistical signiflcance of the latter observation is not very strong. One possibility 
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is that accretion generates a mass quadrupole moment which results in gravitational waves 
and this prevents the star from rotating too quickly [3l[. Upper limits on the generation of 
gravitational waves from this mechanism have already been generated by LIGO 3^ . r-mode 
oscillations can be an alternative explanation for the limiting frequency. The observation of 
these fast-rotating objects provides an important constraint on the extent to which r-mode 
oscillations can spin down neutron stars. Since inferred temperatures of an LMXB's interior 
are much less than for a quark superfiuid, it is pertinent to ask how quark matter inside 
an accreting neutron star affects the r-mode and the star's rotation speed, assuming the 
quark core is not spun out of existence. 

Motivated by these reasons, in this article, we investigate the effect of quark matter on 
r-modes in compact stars. In section [Til we calculate the r-mode frequency for perturbed 
rotating fluids that obey chosen nuclear or quark matter EoS. In section UTTI we review the 
viscosities of neutron matter, ungapped and gapped CFL quark matter, which are inputs in 
a calculation of the damping timescale of the r-mode. In section IIV[ we present results for 
r-mode damping times in strange and hybrid stars, and discuss the critical spin-frequency 
in section |Vl In section Wl\ we remark on observational consequences for such compact 
stars and present our conclusions. 



II. r-MODES IN HADRONIC AND QUARK MATTER 

Various kinds of pulsation modes exist in neutron stars, classified by the nature of the 
restoring force (see Ref. j33| for a review). For example, p-modes are high-frequency (few 
kHz) pressure waves where fluid oscillations are largely radial, while (7-modes are low- 
frequency (~few 100 Hz) density waves driven by gravity which tends to smooth out compo- 
sition and thermal gradients, particularly in proto-neutron stars. The r-mode is intimately 
linked to the rotational properties of the star, and the restoring force here (in rotating stars) 
is the Coriolis force. The modes are termed quasi-normal when they lose energy to gravi- 
tational waves and may be classified by the parity of the fluid displacement vector: polar 
(spheroidal) or axial (toroidal), r-modes are purely toroidal only for the trivial case of a 
non-rotating star where they have zero frequency. In a rotating star, the fluid displacement 
vector corresponding to the r-mode acquires spheroidal components as well, complicating 
the mode analysis. However, it is conventional to assume that the modes in rotating stars 
remain toroidal. It is also conventional to apply the Cowling approximation ^ in which the 
back-reaction of the fluid perturbation on the metric, and hence the gravitational poten- 



tial, is ignored. Previous studies |3J, |35| have shown that including such effects modifies 
the r-mode frequency approximately at the 5% level. Furthermore, if the perturbation is 
assumed to be isentropic, 6P and 6p obey the same EoS as the unperturbed quantities: P 
(the pressure) and p (the density). With these approximations, the r-mode frequency in the 
co-rotating frame, to first order in the rotation frequency of the star Q, is given by js^] 



^ The Cowling approximation is equivalent to neglecting (50o in Eq. ([¥]). We include (50o, however small, so 
we do not use this approximation. 
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Since we are interested in the instability to gravitational wave emission, we restrict ourselves 



to the "classical" r- modes of Papaloizou and Pringle 
observer measures a r-mode frequency of 
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from which it can be deduced that, for / = m > 2, a counter-rotating mode in the rotating 
frame appears as co-rotating with the star to a distant inertial observer. Thus, all r-modes 
with m > 2 are generically unstable to the emission of gravitational radiation and the 
m = 2 r-mode is the first to go unstable. This is the Chandrasekhar-Friedman-Schutz 
(CFS) mechanism js^, 3^- This instability is active as long as its growth-time is much 
shorter than the damping time due to the viscosity of stellar matter. Its effect is to slow 
the rotation rate of a compact star in a short time span of about a year 2^. This possibly 
explains why only slowly rotating pulsars are associated with supernova remnants. The 
r-mode instability might not allow millisecond pulsars to be formed after a supernova or 
an accretion induced collapse of a white dwarf; rather it seems that millisecond pulsars can 
only be formed by the accretion induced spin-up of old, cold neutron stars. However, in 
such cases, the neutron matter inside the neutron star could be in a superfluid state, and 
conclusions for r-mode damping in a superfluid phase are strongly model-dependent |39j. 

Studies of r-mode oscillations provide a unique avenue to sample the star's density profile 
in addition to its mass and radius, r-mode oscillations distinguish between "pure" neutron 
and hybrid or strange stars with the same mass and radius. To leading order in Q, the 
r-mode frequency has no dependence on the EoS, Eq. ([T]). Including second-order rotational 
effects, one finds the following relation between the mode frequency Uj. in the inertial frame 
and n M 
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rot 



niQ = kQ — mQ, 



(3) 



m 



1^2 



irGpo 



where ^2 is obtained explicitly from Eq. (j4]) below, po is the average density of the unper- 
turbed star and G is Newton's constant. Ignoring general relativistic effects, po is related to 
the Kepler frequency as Qk 



|-\/27rGpo |41[. ^2 depends on the density profile of the star, 
as well as density and gravitational perturbations, which are related by the Poisson equation 



for (Newtonian) gravity. We determine K2 from the second-order analysis of r-modes 40 
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where < r < R is the radial coordinate in the star, P22('") is the radial dependence of the 
non-isotropic correction to the density of the rotating configuration to lowest order in Q and 
S^o{r) describes the radial dependence of the change in the gravitational potential due to 
the r-mode oscillation, again to lowest order in Q. The quantity h is the barotropic enthalpy 
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and {dp/dh)o denotes a derivative in the unperturbed star. The function S^Q{r) is deter- 
mined through the solution of the differential equation 
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whose derivation is detailed in Ref. [40|]. The function P22{f) is related to ^22(^)1 the 
change in the gravitational potential due to centrifugal deformation to lowest order in Q, as 
P22(r) = (^)q [$22 (r) - InGpor^], where $22 (r) obeys |4o| 
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enforced via logarithmic derivative matching 
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Along with the condition that the perturbations vanish at r = 0, Eqs. ([6]), ([7]) fix P22{r) 
and S^o{r) throughout the star and subsequently K2 through Eq. (j4]). We determine K2 and 
hence the r-mode frequency for the I = m = 2 mode since this is the lowest / that can couple 
to gravitational waves. We consider four equations of state; the first two are for hadronic 
matter, the next two are for quark matter. 

1. The Akmal-Pandharipande-Ravenhall (APR) hadronic EoS 4^ whose microscopic 
input is based on the Argonne vis nucleon-nucleon interaction j43| which is calibrated to 
deuteron properties and vacuum nucleon-nucleon phase shifts for laboratory energies E'lab 
up to 350 MeV . The inclusion of the Urbana IX three-body force [13] and a relativistic 
boost term 6vh successfully reproduces binding energies of several light nuclei. 

2. To make contact with earlier works, we also employ a polytropic EoS 



Kp 



1+1/n 



(9) 



where i^T is a dimensionful constant for finite values of the polytropic index n. The case 
n = denotes incompressible matter, which softens with increasing n. In Newtonian gravity, 
stable configurations exist only for F = 1 + 1/?^ > 4/3, i.e., for n < 3. We do not include 
general relativistic corrections in the structure or in the analysis of r-modes. Omitting these 
corrections leads to a smaller mass and larger radius for the compact star (see values of 
radius R in Table [IT|l . The effect of general relativistic corrections on the r-mode spectrum 
remains to be analyzed. 

3. An MIT Bag model EoS for charge- neutral self-bound ungapped quark matter, with 
Bag constant B and current quark mass m<j. The Bag constant is chosen such that the 
energy per baryon E/A < 930 MeV at P=0. The pressure is given by 
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For the parameters used in the calculation rriu ~ ma <^ <^ /Xj, we can write in perturba- 
tion 

+ °<'""- (11) 

i=u,d,! 



with quark chemical potential 



fit = fig- QiHe, (12) 



where Qi is the electric charge and /i^ = (/!„ + /id + /is)/3 is the average quark chemical 
potential. Charge neutrality dP/dfie = implies fie ~ ''^1/ i.'^l^q) S^^ 



3^ 3^ 



Ug = Uu + rid + ns 



dP 3/i3 3^^rn2 



The energy density e = —P + Ugfig yields 



471-2 471-2 

Eliminating /ig between Eqs. ( 1131) and ( iMl) gives the EoS 



^ - (14) 



1 , . „x 



- 3.6-4i?)-^v/7^, (15) 

which when compared to the exact EoS, is accurate to about 2% for the entire density range 
of interest in a strange star made of ungapped quark matter. 

4. The CFL phase of quark matter [l| with non-zero quark masses. This phase is 
characterized by a common quark Fermi momenta kp and Cooper pairing among all three 
flavors and colors of quarks. This implies that the number densities of all three quark 
flavors are equal, which enforces electric c harg e neutrality even without electrons 0, 46 
The pressure in the CFL phase is given by [47| 



/ dkk\^Jk^ + m^-^ig) + (16) 



CFL — ^ 



where small contributions from the phonon and pseudo-Nambu-Goldstone bosons have been 
neglected since they have a negligible effect in determining the stellar structure. For realistic 
quark chemical potentials fig, the strange quark's current mass ^ m„ ~ nid can cause 
a mismatch in the quark flavor chemical potentials, but as long as the mismatch is much 
smaller than the superfluid gap A, Cooper pairing with a common Fermi momentum kp is 
still favored. The condition for a common Fermi momenta < 4/igA j48| is satisfled for 
the realistic parameters we use in the calculation. The average quark chemical potential 
f^q = ifJ'u + fJ'd + /Us)/3 determines the total quark number density 

dP 3 

ng = nu + ns + nd = ^ = — (A;|, + 2AVg), (17) 
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which is three times the number density for each flavor. In perturbation where m„ ~ <C 
nT-s <S /Ug, we have 



2 



kF = f^,-^ + 0{ml ml mt), (18) 



with pressure and energy density 

,4 Q,,2 



Pc..-^-5 + ^(4A^-mD, (19) 

Just as for ungapped quark matter, at leading order in the quark mass expansion, we have 

PcFL ^ -{e-AB) + ^t^±-^VT^. (20) 
o on 

This approximate form based on the expansion in m^/fig works better than expected due to 
the accidental combination (4A^ — m^)//Xg which is numerically small. 0{m'^) corrections 
simply lead to a renormalization of the bag constant [lo| and can be ignored. Indeed a 
comparison to the exact EoS verifles that the approximate form Eq ( |20l) is accurate to 
within 0.1% for the entire density range of interest in a strange star. 

For hybrid stars, we employ the APR EoS for the hadronic phase and either the MIT 
Bag or CFL EoS. Without modiflcations, a combination of the equations of state above 
does not generate a stable stellar conflguration which contains a signiflcant amount of quark 
matter and which also gives a sufficiently large maximum mass when general relativistic 
effects are included. For this reason, we make two additional modiflcations. Firstly, we 
compute the hybrid star proflles using general relativity, even though the expressions for 
r-mode properties are derived in the Newtonian limit. While a more consistent approach 
is desirable, the present approach is not unreasonable because it is the composition of the 
star, not the gross structural details, which determines the conclusions about the stability of 
the r-mode. Secondly, we add a term — 3c/i^/(47r^) to the pressure, which characterizes the 
0{as) perturbative QCD correction to the pressure of a non-interacting quark gas [lo|, 
with < c < 0.3. 

For hybrid stars, the bag constant and the parameter c are chosen such that transition 
density between hadronic and quark matter comes out to be about twice the nuclear satu- 
ration density and the maximum mass in general relativity is 1.7 Mq. For ungapped quark 
matter, B = 75 MeV/fm^ and c = 0.083, and for gapped quark matter, B = 106 MeV/fm^ 
and c = 0.12. ^ Unless strange quark matter is absolutely stable at zero pressure, lower tran- 
sition densities are incommensurate with the nonobservation of deconflned quarks in nuclei. 
Thus, the procedure outlined above effectively maximizes the amount of quark matter for 
a given hadronic equation of state in a hybrid star. Conflgurations with higher transition 
densities will have smaller quark matter cores, with results that lie between those for pure 



^ CFL matter has larger pressure at a fixed density than ungapped matter because of the gap; consequently, 
the Bag constant value must increase to ensure that the phase transition occurs at the same density in 
both cases. 
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hadronic stars and hybrid stars with maximized quark content. Fig. [T] displays the typical 
density profile for hybrid stars containing ungapped or CFL quark matter for stars with 
mass 1.4Mq. 

We assume that the surface tension between the quark and hadronic phases is large 
enough so that only one surface is created. An estimate of the critical surface tension 
and the microscopic nature of the minimal CFL- nuclear interface is given in Ref. 47|. The 
assumption of a sharp boundary is made in the interest of simplicity, in which case a Maxwell 
phase construction is sufficient to describe the surface, with pressures and chemical potentials 
being continuous across the surface but density being discontinuous. A density discontinuity 
is clearly an idealization of a microscopically thin surface region where, for example, the 
number density will continuously transition from the larger value in the quark phase to the 
smaller value in the hadronic phase. This procedure is a good approximation if the surface 
tension is large enough that the energy gain provided by separating matter into a mixed 
phase is smaller than the surface energy cost. 
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FIG. 1: The density profile of hybrid stars containing ungapped (APR+Bag) and gapped 
(APR+CFL) quark matter that are considered in this work. The transition densities were fixed at 
0.32 fm~^. The vertical lines demarcate the homogeneous quark and hadronic phases. 



For the equations of state discussed above, results for the r-mode frequencies are listed 
in Table 1. We confirm earlier findings for r-mode frequencies in neutron matter described 
by polytropic EoS, viz., that the r-mode frequency for polytropes increases with increasing 
softness of the EoS (higher polytropic index n). For gapped or ungapped quark matter, the 
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r-mode frequency is similar to matter with a polytropic index ^ 1.6. The r-mode frequency 
is robust to variations in the gap < A/(MeV) < 100, as long as other parameters are 
adjusted such that quark matter remains self-bound. 

TABLE I: H2 from Eq. Q and m=2 r-mode frequency from Eq. ^ for various equations of state. 
We choose the central density such that the star's mass is 1.4Mq in each case (unless otherwise 
specified) while the rotation frequency is chosen to be = Qk, where is EoS dependent. 
For quark matter in strange stars, we choose B =80 MeV/fm^, 771^=100 MeV and A = MeV 
(ungapped) or A =100 MeV (gapped). For hybrid stars, we choose B = 75 MeV/fm ^ and 
c = 0.083 (ungapped) and B = 106 MeV/fm^ and c = 0.12 (gapped). The fractional correction 
over the 0*'^-order EoS-independent result loI^^ given by Eq. ([2]) is also displayed. 



EoS 


K2 


LOj. ^(ms) 


\uJr/^\ 




n = 1 Poly trope 


0.298 


0.155 


1.215 


0.088 


n = 2 Polytrope 


0.148 


0.147 


1.275 


0.044 


APR 


0.316 


0.180 


1.208 


0.094 


APR (M=1.8 Mq) 


0.266 


0.161 


1.228 


0.079 


Bag 


0.198 


0.106 


1.255 


0.059 


CFL 


0.195 


0.099 


1.256 


0.058 


APR Bag (A/=1.7 Mq) 


0.083 


0.097 


1.301 


0.025 


APR CFL (M=1.7 Mq) 


0.099 


0.100 


1.294 


0.029 



Since the r-mode is really a quasi-normal mode, the mode frequency will generally acquire 
an imaginary part on account of the energy transmitted in gravitational waves and through 
dissipative forces in the fluid. To assess the r-mode instability and critical rotation frequency 
of compact stars, a comparison of the viscous damping timescale with the mode-growth 
timescale due to gravitational wave emission is required. We proceed to a discussion of the 
role of viscosity in r-mode damping. 



III. VISCOSITY OF HADRONIC AND QUARK MATTER 

The temperature-dependent bulk viscosity ( and shear viscosity r] of the fluid suppress 
r-mode growth. For hadronic matter, we use simple power-law fits derived in Ref. [Hi which 



were shown to faithfully reproduce results of microscopic calculations of the viscosity [51 
For non-superfluid neutron matter, the dominant contribution to the bulk viscosity comes 
from the modified Urea process: n + n^n + p + e~ + p, 



-2 



g/(cm s), (21) 

where the dimensionless quantities pis = pj (10^^ g/cm"^) and Tg = T j (10^ K). The dominant 
contribution to shear viscosity comes from nn scattering and is given by 

77„„ = 2 X lO^V^^'Tg"' g/(cm s). (22) 
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At low temperatures neutron matter is likely to be in a superfluid state. At higher den- 



sities in the stellar interior, neutrons pair in the state with ~ 10^ K j52| while they 
5*0 state with Tc ~ 5 x 10^ K jss], 13] in the outer layers. At temperatures be- 



pair m a 

low the critical temperature Tr for superfluidity, the viscosity contributions from modifled 
Urea processes are suppressed |55|] so that the phonons and Nambu-Goldstone modes corre- 
sponding to spontaneous breaking of rotational invariance ( "angulons" ) in the phase 5^ 
dominate the viscosity. Ideally, the viscosity estimates [Eqs. f l21l) . fl22p ] for neutron matter 
at T <^ Tc should incorporate the effects of superfluidity. We neglect these effects in the 
interest of simplicity, since the focus is on the quark phases. We note, however, that Lind- 
blom & Mendell |39|] have shown that K2 (and hence the r-mode frequency) for a superfluid 
neutron star described by the APR EoS differs by less than 0.1% from the non-superfluid 
case. Furthermore, in computing the viscous damping timescales, these authors have used 
electron scattering off magnetized vortices, termed mutual friction 57|], as the dominant 
source of bulk viscosity, and electron-electron scattering as the dominant source of shear 
viscosity. The results depend strongly on the model of superfluidity in the core and intro- 
duce additional complications in the analysis. We therefore postpone the practical issue of 
r-modes in superfluid neutron matter to future work. 

For ungapped quark matter, the bulk viscosity is determined mainly by the weak process 
d + s ^ u + s jll], while leptonic contributions can become important at high temperature. ^ 
Since we are only examining stability and not attempting to describe the full evolution of 
the r-mode, we use the approximate expression for the bulk viscosity [Eq. (16) of Ref. (58| ] 
which is appropriate for small oscillations of the fluid and when 27rT = /x^ — /^d (this 

holds in the range of temperatures of interest in this work). We have 



l-[l-exp(-v^TV)] 



( Ji) G'^sin^^.cos^^.mt/ilT^ 



6.66 X 10 



20 



_/id_y 

MeV/ 



N 4 



g/ (cm s^ 



3.57X 10-Mr^) ( 1 + 



MeV/ 



(23) 



where u = Ur = 2tt/t is the frequency of the perturbation, Gp the Fermi constant and 
9c the Cabbibo angle. In the bulk viscosity expression, the second term inside the big 
square bracket is numerically negligible for the parameters of the calculation. Further, this 
"transient" term drops out when averaged over oscillation cycles leaving only the overall 
factor in front, similar to the CFL bulk viscosity expression in Eq. (!25l) . 

The shear viscosity of ungapped quark matter is dominated by quark-quark scattering in 



At high temperatures T > 10 K, the semi-leptonic processes d(ov s) —f u + e~ + and u + e~ —>■ 
d{oT s) + Ve can also contribute significantly to the bulk viscosity 59|. However, for our parameter choices 
of -nis = 100 MeV and Wr=l kHz, this contribution is negligible. 
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QCD. Neglecting small QED effects, we have 6^ 



/ (cm s). 



(24) 



In the superffuid CFL phase, the quarks are gapped. The massless phonons and thermally 
excited light pseudo-Nambu-Goldstone bosons determine the thermodynamic and hydrody- 
namic properties. Unlike the mesons of the QCD vacuum, the mass ordering of the CFL 
pseudo-Nambu-Goldstone bosons is "reversed" and the neutral kaon is the lightest me- 
son f6l| . The contribution to bulk viscosity from phonons alone has been calculated in 



Ref. 62. However, the dominant contribution to bulk viscosity comes from flavor chang- 
ing K'^ decay, viz., via weak equilibrium processes <-> 00 and (pK^ involving the 

It should be noted that in contradistinction to normal matter, 



massless phonon [63|, 164 



superfiuid hydrodynamics is described by more than one bulk viscosity. Additional viscosi- 
ties, commonly labeled Ci and (^3, enter the dissipative terms in the stress-energy tensor, 
and they multiply factors of Vs — Vn, the difference in the speed of the "superfiuid" and 



"normal" component in the two- fluid model of superfluidity (65|, l66|. The superfiuid com- 



ponent is non-viscous while the normal component consisting of the phonons and other low 
energy excitations such as pseudo-Nambu-Goldstone modes is viscous. Due to the large 



thermal conductivity of neutron and CFL matter [1^, |67j, we may approximate the star 
to be isothermal. Consequently, the second sound associated with a non-zero Vg — Vn and 
temperature gradient would be damped, and contributions from (i and (3 can be ignored 
in the leading order of this approximation. For the calculations, we employ only the bulk 
viscosity G (which becomes the usual bulk viscosity ( in normal matter) from neutral kaon 
decay B: 



Cro = C 



c 



2 ' 



UJ' + 7efr 



(25) 



"-It-", if)' 



7efr = iK 



dflq 



2 ' 



duq driK \ 



where as before u = Uj. = 2-k /t^ and 7efr is the effective kaon width [6J]. The number den- 
sities for quarks Uq and neutral kaons uk are determined from the thermodynamic pressure 
P = Pcfl + Pr where the contribution from kaons was added to the pressure Pcfl from 
Eq. (fT6|) . We choose parameters such that the effective kaon chemical = (m^— m^)/(2/ig) 
is smaller than the kaon mass iriK ensuring that there is no kaon condensation. To com- 
pute the bulk viscosity, one considers small departures from equilibrium characterized by a 
non-zero chemical potential 5^k and writes 



Pk 



27r2 



dkk"^ In 



1 — exp 



E 



R 



5ji 



R 



T 



(26) 
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The barred quantities and derivatives in Eq. (125!) correspond to equilibrium values S^k = 0. 
The width 7eff ~ 'Jk is determined from the kaon decay rate [g^] 

f driK \ rforward(5/iir = 0) , . 



forward ~ /- off2 r _ , / /tt^x . i 2 ' 



gySTT \^ ^-f) [expp/(v^T)-l]' 

2 eff2 



cff 



where y^d; are the usual CKM matrix elements. Perturbative estimates for the decay 
constants yield [gH 



^2 21 - 8 In 2 /i; 



2 



1^2^' ^^^^ 



f2 



3 /ia 



4 27r2' 



In the bulk viscosity calculation we use A = 100 MeV, = 5 MeV, = 7 MeV, 
ms = 100 MeV and B = 80 MeV/fm^ to construct a strange star in the pure CFL phase 
with mass 1.4 Mq. A self-consistent determination from Eq. ( l20l) gives fiq ~ 300 MeV 
and /i^ ~ 17 MeV. Unfortunately, at these moderate densities the kaon mass rriK is not 



accurately determined but it is expected to be about 15 — 25 MeV [68|]. If there is a region 
in the star where Sm = mx — < 0, then modifications to the equation of state and 
viscosities due to kaon condensation would be required j69|, but to study a concrete case, 
we treat rriK as a free parameter and choose 6m = rriK — ^j^^k — ^ MeV, avoiding the issue 
of kaon condensation. The CFL bulk viscosity drops rapidly with decreasing T at higher 
values of 5m. 

The shear viscosity in the CFL phase has been calculated from phonon scattering 00 ^ 



00 |70|. The contribution from the medium modified photon was shown to be small for 
temperatures where the phonon viscosity is of interest [70|, therefore it is not included in 
the present calculation. A variational estimate of the linearized Boltzmann equation relevant 



for shear viscosity gives [70|] 

r/^ = 3.745 X 10^ (^)'^9"' g/(cm s). (29) 

The shear viscosity falls steeply with increasing temperature and consequently has a small ef- 
fect on the viscous damping of the r-mode at high temperature, see Fig.[2]and Eqs. (!36|) . (!39ll . 
Although the shear viscosity is undoubtedly more important at low temperatures, the con- 
tribution from the phonon becomes irrelevant at very low temperatures. This is because 
being a Nambu-Goldstone mode, the phonon is derivatively coupled where for typical in- 
teraction terms idocf) = E^(f) ~ T0 for thermally exited phonons. At very low temperatures 
T <^ fig, then, the phonon scattering cross section is suppressed and the mean free path 
increases. The hydrodynamic expression for the shear viscosity is only relevant when A,^ is 
much smaller than the star radius R, otherwise the phonon travels through the star without 
collisions responsible for shearing the fluid flow. An order of magnitude estimate of the rel- 
evant phonon mean free path can be obtained from the kinetic theory relation t]^ ~ 
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10"^ 10'^ 10 ' „ 1 10 10^ 10^^ 10"^ 10 ' „ 1 10 10^ 
T, = T/IO' K T, = T/IO' K 



FIG. 2: The temperature dependence of the bulk and shear viscosity of neutron matter from 
Eqs. ([2T]) . (p2]) (left panel) and that of ungapped and gapped quark matter from Eqs. ([23j) . (p^ 
and Eqs. (j25p . (j29p respectively (right panel). For neutron matter, the energy density is chosen to 
be 1.5 X 10~^ MQ/km'^ and kQ, is fixed at 1000/s. For quark matter, the quark chemical potential 
is chosen to be 310 MeV, =100 MeV with A=100 MeV for CFL matter, and r = In/uJr is fixed 
at 0.001 s. 



where is the phonon number density and p ~ 2.7T/v is the thermally averaged phonon 
momentum with phonon speed v = 1/ -\/3- This length scale quantify the distance over 
which the photon must travel to generate shear viscosity due to collision. We estimate 

A . ^ ^ 2.1 X 10^ ( km. (30) 

^ pn^ V 300 MeV/ ^ ^ ^ 

For typical quark chemical potential fig ~ 300 MeV, the phonon mean free path will 
exceed the stellar radius at temperatures T < 10^'^ K. The calculated damping times for 
the CFL phase are valid only above this temperature. For T < lO^'' K, contributions other 
than from phonons are expected to be relevant for shear viscosity but these have yet to 
be estimated. For this preliminary work on the CFL r-modes, we have included only the 
phonon contribution. 

For hybrid stars, we ignore possible bulk viscosity due to the quark-hadron interface j7l| , 
and simply utilize the viscosity expressions above for each phase. This approximation implies 
that the viscosities for hybrid stars may be underestimated compared to a realistic treatment 
of the interface/crust. 



IV. VISCOUS DAMPING OF r-MODES 



The energy of the r-mode is dissipated according to 7^ 



dE 



(31) 
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where 5a = VaSv"" is the volume expansion due to the r-mode and Saab = |(Va5t'fe + Vb^fa — 
l^afoVc^t'^) is the shear tensor. The first term is the energy radiated in gravitational waves 
to lowest order in Q with 6Jmm being the current multipole 



R 



(m.+l) 



(32) 



For m > 2, uJrot < rriQ, so that the r-mode energy grows with gravitational wave emission, 
triggering the instability. The timescale r associated with growth or dissipation is given by 



1 



1 

2E 



The energy E of the r-mode is given by 

E=—(m + lf (2m + iV.R^n'^ / drp ( -) 
2m Jq \RJ 



2m+2 



(33) 



(34) 



where we have dropped a proportionality constant that determines the amplitude of the 
r-mode, since it cancels in the evaluation of the damping timescale. Explicitly, the gravita- 
tional radiation time scale is \71 



tgw c2™+3 [(2^+1)!! 

while the shear viscosity time scale is [tI] 

1 {m — l){2m 



m + 2 



m 



2m+2 



drpr 



2m+2 



(35) 



lo drpr 



2m+2 



drrjr 



2m 



(36) 



We compute the bulk viscosity timescales using a simplified expression for the volume 
expansion, which is based on approximating the Lagrangian perturbation of the fluid by an 



Eulerian one 72 



6 



a 



Sp 

p ' 



(37) 



As shown in Ref . [40| , while this approximation is only good within an order of magnitude for 
the bulk viscosity damping timescale, the critical angular velocity, which can be constrained 
by observations, is hardly modified because the bulk viscosity has a steep temperature 
dependence. Thus, the critical frequency curves presented in this work are accurate, even 
though the bulk viscosity timescales are only approximate. 

Carrying out angular integrations for the bulk viscosity part of Eq. flM]) and using 
Eq. ([ST]), one obtains 




^ =-2 



(2m + 1)! 
2m + 3 



dr{27rr^)C 



dp 

(1 + m)2p2 



r \ m + l 



-| 2 



(38) 



where po{r) is the density of the unperturbed star. The bulk viscosity damping timescale 
follows from Eqs. (I33l). (IMj) and (|38|). 
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For a fixed temperature T=10^ K, and for stars rotating at tlie Kepler frequency Q = fix, 
all three damping timescales are displayed in Table [Tll The timescales are also plotted as 
a function of temperature in Fig [3l We observe from Fig. [3] that in non-superfiuid neutron 
matter, shear viscosity damps the r-mode rapidly only at very low temperatures T < 10^ K, 
while bulk viscosity damps it only at very high temperatures T > 10^° K. In normal quark 
matter, the damping timescale due to bulk viscosity has a minimum around T ~ 10^ K, 
corresponding to a maximum in the bulk viscosity. This behavior leads to a high-temperature 
stability window. In CFL quark matter, the damping timescale due to shear viscosity at 
low temperature is large. However, the phonon contribution considered here is not reliable 
below about 10^° K as mentioned earlier. Other contributions, such as from light-by-light 



scattering [70|, could become important in this regime and damp out the r-mode. Keeping 
this caveat in mind, and noting that the damping timescale due to bulk viscosity is still 
large except at high temperatures T > 10^^ K, the r-mode in CFL matter is unstable only 
in a small window around high temperatures T ~ 5 x 10^° K, in marked contrast to either 
ungapped quark matter or normal neutron matter. It would be interesting to see how the 
width of this window changes upon inclusion of relevant shear viscosities at temperatures 
T < 10^° K. 




FIG. 3: The temperature dependence of the bulk and shear viscosity damping timescales for 
neutron matter (left panel) and for strange and hybrid star matter (right panel) . The rotation rate 
is the Kepler frequency i^x- 



V. CRITICAL ROTATION FREQUENCIES 



The critical rotation frequency = 1/tc of neutron/strange/hybrid stars can be deter- 
mined by the criterion that at this frequency, the fraction of energy dissipated per unit time 
exactly cancels the fraction of energy fed into the r-mode by gravitational wave emission: 



1 



' 1 



1 



1 



0. 



(39) 
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TABLE II: Damping times for various equations of state. The central energy density pc in each 
case is chosen to yield a 1.4 Mq star (unless otherwise explicitly stated) and the temperature 
is fixed at T=10^K. Tf is the critical or cumulative r-mode damping timescale given by Tf = 
(1/r^ + l/r,; + 1/tgw)~^. The star is taken to be rotating at the Kepler frequency VIk- Energy 
densities are given in units of 10"^ Mo/km^. 



EoS 


Pc 


i2(km) Q 


;^(kHz) 


r^(s) 


Tr;(s) 


^Gw(s) 


rf(s) 


n = 1 Polytrope 


5.63 


12.50 


5.31 


2.47 xlO^° 


2.48 xlO^ 


-5.22 xlO^ 


-5.22 xlO^ 


n = 2 Polytrope 


19.5 


12.50 


5.31 


6.20 xlO^° 


7.28 xlO'^ 


-1.29 xlO^ 


-1.29 xlO^ 


APR 


3.60 


13.74 


4.61 


2.48 xlO^ 


3.25 xlO^ 


-9.54 xlO^ 


-9.54 xlQi 


APR (M = 1.8 Mq) 


3.97 


14.06 


5.05 


1.79 xlO^ 


3.03 xlO*^ 


-3.58 xlO^ 


-3.58 xlO^ 


Bag 


4.30 


9.89 


7.54 


7.33 xlQi 


2.78 xlO^ 


-8.28 xlO° 


-9.34 xlO° 


CFL 


3.23 


10.71 


6.70 


1.74 xlO^o 


5.38 xlO-i 


-7.48 xlO"^ 


1.92 xlO° 


APR+Bag (M=1.7 Mq] 


13.5 


10.25 


7.93 


1.04 xlO^ 


2.38 xlO^ 


-7.96 xlO° 


-8.62 xlO° 


APR+CFL (M=1.7 M0 


) 10.6 


10.40 


7.69 


1.72 xlO^ 


1.29 xlO"^ 


-8.66 xlO° 


-8.67 xlO° 




T (K) T (K) 

FIG. 4: The critical frequency in units of as a function of temperature for neutron stars, 
strange stars and hybrid stars. For quark matter, we choose ms=100 MeV and A=100 MeV. The 
bag constant i?=80 MeV/fm^ for strange stars, while B=11Q (150) MeV/fm^ for hybrid stars with 
ungapped (CFL) quark matter. The box represents typical temperatures (2 x 10^-3 x 10® K) and 
rotation frequencies (300-700 Hz) of the majority of observed LMXBs assuming = 5500 Hz. 



Stable rotation frequencies at any temperature will be limited by the smaller of the 
critical frequency or the Kepler limit f2/f2/f=l. As depicted in Fig. |U the region above the 
temperature-dependent VLc curve is unstable to r-mode oscillations and the star, if it enters 
this region, will be spun down rapidly to < VLc- 

Normal neutron stars are generally unstable to r-modes except at very high or very low 
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temperatures, where the bulk and shear viscosities respectively are effective at damping 
r-mode oscillations. The width of the instability window is not very sensitive to stellar 
mass, but is noticeably smaller for a softer equation of state (larger polytropic index). The 
box marked "LMXBs" is representative of typical LMXB core temperatures [73| 2 x 10^ 
to 3 X 10^ K, with rotation frequencies between 300 and 700 Hz. Strange stars made of 
ungapped quark matter are more stable than their neutron star counterparts in a window of 
temperatures 10* K to 5 x 10^ K where the bulk viscosity damping timescale in quark matter 
is quite small. Strange stars composed of CFL matter are unstable to r-mode oscillations 
in the range lO^'' K< T < 10^^ K. Hybrid stars with ungapped quark matter are much 
like ungapped strange stars, except that the bulk viscosity from the outer shell of hadrons 
damps the r-mode oscillations at higher temperatures, r-mode oscillations in hybrid stars 
with gapped quark matter are unstable in the range 10^ K < T < 10^° K. The results for 
hybrid stars are much like the corresponding strange star model because we have maximized 
the amount of quark matter by fixing the deconfinement transition density to be small. If 
this transition density is larger, the critical curves for hybrid stars look more similar to those 
in the left panel of Fig. m 

VI. CONCLUSIONS 

We have calculated the frequency and viscous damping timescale of the m = 2 r-mode 
for compact stars made of hadronic matter, quark matter, and for hybrid stars. To our 
knowledge, this is the first attempt at studying r-mode characteristics for pure or hybrid 
stars with color superconducting quark matter, particularly the CFL phase. In absolute 
terms, the r-mode frequency is about 50% larger for strange stars in comparison to neutron 
stars of the same mass with a polytropic EoS, assuming both are rotating at their respective 
Kepler frequency. The ratio of the r-mode frequency to the star's rotation frequency is 
larger for strange stars and hybrid stars; this is due to the softer equation of state. It is also 
evident that EoS effects on the mode frequency are smaller for softer equations of state. 

The dominant bulk and shear viscosities of the normal hadronic phase, the ungapped 
quark phase as well as the gapped CFL quark phase, which was derived recently, are used 
to estimate viscous damping timescales for homogeneous and hybrid stars. The current 
results confirm established results for neutron stars; viz. that since neutron stars are stable 
against the r-mode instability only at very high (T > 10^" K) or very low temperatures 
(T < 10^ K), they would be spun down rapidly by the r-mode instability shortly after their 
birth at MeV temperatures, when large neutrino losses cause rapid cooling. The predicted 
critical frequency for LMXBs with hadronic neutron stars is close to or inside the region 
given by observations and LMXBs. A more careful analysis taking the effect of superfluidity 
in neutron matter into account would provide a clearer picture on the presence of r-modes 
assuming no quark matter is present. 

Strange stars with non-superfluid quark matter display a stability window between 
10* K < T < 5 X 10^ K where they can spin at a substantial fraction of the Kepler frequency. 
If LMXBs contain strange stars, then on statistical grounds, we should expect to observe 
some of them spinning close to their Kepler frequency, unless some mechanism not connected 
to r-modes is responsible for limiting their frequency or the quark phase is more complicated 
than considered here (ungapped or CFL). If newly-born superconducting strange stars occur 
in nature, they may still spin down rapidly at birth, provided the cooling timescale is slower 
than the spin-down timescale. A consistent analysis of spin-down and cooling is currently 
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under investigation and will be reported in separate work. Hybrid stars tend to be more 
resilient against the r-mode instability, except for a temperature range which is dependent 
on the quark EoS. This effect is not surprising: in hybrid stars, r- modes which would nor- 
mally be undamped in hadronic matter are damped by the quark core and vice versa. We 
conclude that hybrid proto-neutron stars are not likely to be damped by r-modes at all. 
This work also has important implications for LMXBs if they are hybrid stars containing a 
significant amount of deconfined quark matter in the ungapped or CFL phase, or if LMXBs 
are strange stars with a thin accreted hadronic crust. Based on the estimates of viscous 
damping timescales, the presence of quark matter in LMXBs implies that either they must 
contain rapidly rotating stars that have not yet been observed, or have a more complicated 
phase structure. A similar explanation regarding the LMXB rotating at 1122 Hz was put 
forth recently [t^ ]. 

To make the results more concrete, future work will have to address several assumptions 
more closely. We assumed that the modified Urea process provides the dominant contribu- 
tion to bulk viscosity while neutron-neutron scattering determines the shear viscosity. In 
practice, in the core of neutron stars, for T < ~ 10^ K, neutron matter is likely to be in 
a triplet-paired superfluid state. This will suppress the modified Urea rate [HHl, rendering it 
insignificant at T -C Tc. Singlet pairing among neutrons and/or protons has a larger critical 
temperature Tc ~ 5 x 10^ K jSJ] and can occur in the crust with similar consequences. Under 
these conditions, it is more realistic to include alternate sources of viscosity along the lines 
of Ref. 13] • It would also be useful to explicitly check the contributions from the additional 
viscosities Ci and ^3 that are inherent in a 2-fluid model of superfluidity. 

We also assume that in the CFL phase, the shear viscosity is from phonons alone, whereas 
below T ~ 10^° K, the damping timescales are not reliable as the phonon mean free path 
exceeds the stellar size. At these low temperatures, other contributions such as from photons 
would be more relevant and need to be explicitly computed and included for an accurate 
treatment of r-modes in CFL matter. 

In the current work, we considered a simplified model for a neutron/strange star in this 
work, ignoring realistic features such as a crustal layer, rubbing friction at the quark-hadron 
interface and possibly the role of other contributions to the bulk and shear viscosity at very 
high or very low temperatures. Neutron star crusts composed of hadronic matter typically 
damp r-mode oscillations further. Early work in Ref. [t^ suggested that damping timescales 
were decreased by as much as a factor of 10^ — 10^ because of the additional viscosity provided 
by the neutron star crust. This conclusion has been softened in some later works: Ref. [t^ 
concluded that the associated decrease in the damping timescale was a factor of 10^ — 10^ 
smaller, and Ref. [TTj has pointed out that the damping timescale is sensitive to the crustal 
composition. The final result may well be that r-modes are not as important a mechanism 



as believed for spinning down fast-rotating compact objects. Ref. [78[ has also stressed the 
importance of updating the multi-component nature of equation of state in order to fully 
assess the effects of superfluidity. However, this does not change the central conclusion, viz., 
if r-modes do play a role in spinning down neutron stars, the present work demonstrates 
that this effect is much less efficient in neutron stars containing ungapped or completely 
gapped quark matter. 
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